Virtues and limitations of the truncated 
Holstein-Primakoff description of quantum rotors 



Jorge G. Hirsch, Octavio Castanos, Ramon Lopez-Pena, and 
Eduardo Nahmad-Achar 

Institute) de Ciencias Nucleares, Universidad Nacional Autonoma de Mexico 
Apdo. Postal 70-543, Mexico D. F., CP. 04510 

E-mail: hirschOnucleares . unam.mx 

Abstract. A Hamiltonian describing the collective behaviour of N interacting spins 
can be mapped to a bosonic one employing the Holstein-Primakoff realisation, at the 
expense of having an infinite scries in powers of the boson creation and annihilation 
operators. Truncating this series up to quadratic terms allows for the obtention of 
analytic solutions through a Bogoliubov transformation, which becomes exact in the 
limit N — > oo. The Hamiltonian exhibits a phase transition from single spin excitations 
to a collective mode. In a vicinity of this phase transition the truncated solutions 
predict the existence of singularities for finite number of spins, which are spurious. 
However, renormalisation allows to extract from these divergences the exact behaviour 
of relevant observables with the number of spins around the phase transition, and 
relate it with the class of universality to which the model belongs. In the present work 
a detailed analysis of these aspects is presented for the Lipkin model. 
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1. Introduction 

The dynamics of N two-level systems is well described by the Lipkin Hamiltonian 
PP. Born in nuclear physics, it has found extensive use in quantum optics, in the 
generation of squeezed states [2], multipartite entanglement [3], two-mode Bose-Einstein 
condensates [3], and monomolecular magnets [5]. It represents an approximation to 
ferromagnetic Ising models [6] , exhibiting a second-order phase transition in the limit of 
large number of particles which is well described by mean field techniques pj El E]- In 
most Holstein-Primakoff realisation is employed which, when truncated, provides 

analytical solutions in the thermodynamic limit [TOj [IT] . For a finite number of atoms, 
observables like the ground state energy, the energy gap and the number of excited 
atoms exhibit a singular behaviour at the phase transition, going to zero or to infinity, 
while numerical calculations show that they should remain finite. 

Similar singularities arise when using Holstein-Primakkof description of the rotor in 
more complex systems, as for example in the Dicke model [12] , where N two level atoms 
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are coupled with one electromagnetic mode in a cavity JT3l [TH [15] . While the spurious 
character of these singularities has been exposed previously [161 E]> it seems necessary 
to discuss in detail the situations in which the prediction of the truncated Hamiltonian 
cannot be employed to describe finite systems, and also how useful information can be 
extracted from the spurious divergences through the use of renormalisation techniques 

In this work the mean field description of the Lipkin model is obtained using 
the Holstein-Primakkof realisation of the quantum rotor. It is shown that coherent 
Heisenberg- Weyl states [18] provide a mean field description of the system, exhibit the 
existence of different phases in the parameter space, associated with different qualitative 
descriptions, and phases transitions between them. Analytical expressions are obtained 
for the ground state energy and for the expectation values of the number of photons and 
of excited atoms, which are intensive and are shown to represent the thermodynamic 
limit. A truncated Hamiltonian is built which has analytic solutions, and provides 
the eigenstates of the system far from the phase transitions. Close to them, the 
observables exhibit the singular behaviour mentioned above. The spurious character 
of these singularities, and the way in which relevant critical exponents can be extracted 
from them are discussed in detail. 

2. The Lipkin Hamiltonian 

The Lipkin Hamiltonian [1] describes the collective behaviour of N spins or two-level 
atoms, with energy separation e, which interact by scattering pairs of particles between 
the two levels. In the quasi-spin formalism it has the form 

H = eJ z + + pi (1) 

where J x , J y , J z are the three components of the angular momentum operator, with the 
usual commutation relations [Jj, Jk] = itjkiJh an d 7x>7y are the coupling strenghts. 

2.1. Mean Field description 

In the literature it is customary to employ the Holstein-Primakoff representation of the 
angular momentum operators [19| 

J + = VN#fI§ , J_ = rtf^t, J. = (2) 

where J + = J x + iJ y , J_ = J x — iJ y , the Bose operators b\b obey the commutation 
relation [b,b*] = 1, and the vacuum of the bosons satisfies b\ ) = 0. Making these 
substitutions into H, Eq. Q, the bosonic Hamiltonian is built. Care must taken in 
using the commutation relations to move the creation operators (outside the square 
roots) to the right, and the annihilation ones to the left. 

The mean field description of this Hamiltonian is easily obtained employing as a 
trial state the Heinseberg-Weyl coherent state \a), which is the eigenstate of the bosonic 
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annihilation operator b\a) = a \a>), where a = pe 1 ^ is a complex numbers, with p > 
and < < 2n. 

The expectation value of the Hamiltonian for this co herent s tate provides the energy 
surface. Employing the approximation sjl — ^|/3) ~ y 1 — [20], which becomes 

exact in the thermodynamic limit, when the number of atoms goes to infinity, it takes 
the simple form [2T] 



+ 



2 




(«\H\n)=e[? _^-)+ 2,r) (3) 

7x - 



For a given set of Hamiltonian parameters 6,7^,7^, the values p c , C which minimise 
this expression provide the mean field wave function. They are obtained by solving the 
equations for the energy surface critical points 
d(a\H\a) = Q d(a\H\a) = Q 
dp ' <90 

The solutions of these equations associated with the minima of the energy surface are 

Pc = 0, 0c undetermined, if 7^ > 7 C and 7 y > 7 C , 

(l - J ) , 0c = 0, 7T, if 7^ < 7c and ^ x < 7^ (5) 
Pc = Jf (l - ^) , 0c = 0, 7T, if 7 y < 7c and 7^ > 7 y 

with 7c = — e. The first case, 7^ > 7 C and 7 y > 7 C , defines the normal region (I), where 
in the ground state all atoms are in their lowest energy state. In the first deformed 
region (II), where 7^ < 7 C and 7 X . < j y , it is energetically favoured to collectively excite 
all atoms, and the ground state is doubly degenerate, as there are two critical phases 
0c = 0, Ti. The second deformed region (III), where 7 y < 7 C and 7 X > j y , is symmetric 
to the first one: it is obtained by interchanging J x •f-)- J y . In what follows we will restrict 
our analysis to regions I and II. 

Employing these critical values, explicit analytical expressions can be found for 
intensive observables: the ground state energy per atom e gs = (a c \H\a c ) /N, and the 
fraction of excited atoms n e = 2p 2 c /N: 

7c Jx + l y 
2 47V 

7 c + ll , (7x + 7c)(7y + 7x) = j _ 7c 

4 7 , 87V 7x ' e 7, 

Although they were obtained using different techniques, these mean field 
expressions, shown as dashed red lines in Fig. [TJ exactly coincide with those presented 
in [31 EJ- They reproduce reasonably well the numerical results even for a relatively 
small number of particles, as can be seen in Fig. [TJ where the numerical diagonalization 
of the Hamiltonian matrix H is shown as a continuous green line, for N= 10 and 40 
atoms, for 7 C = — 1 and 7 y = 1. At this point it is worth to emphasize that the mean 



£ gs = 77 H TWr~i n e = 0, region I 

£ gs = -^rr^ + -■> n e = 1 - — . re g lon IL ( 6 ) 



The truncated Holstein-Primakoff description of quantum rotors 



4 




Figure 1. (Color online) Left: Ground state energy per atom, calculated through 
numerical diagonalization of the matrix Hamiltonian H for N= 10 (continuous green 
line), and from Eq. ^ (dashed red line). Right: Fraction of excited atoms calculated 
through numerical diagonalization of the matrix Hamiltonian H for N= 40 (continuous 
green line) and from Eq. ^ (dashed red line). Both plots were calculated using 
7 C = —1 and j y = 1. Peaked blue lines show the spurious results obtained with the 
truncated Hamiltonian (see next section). 



field description becomes exact in the thermodynamic limit, where only the first terms 
in the energy per atom remain finite. The peaks in the continuous blue lines at ^ x = — 1, 
which goes to infinity for the fraction of excited atoms, is similar to the one shown in 
Fig.l of Ref. [T33 - It is a spurious result coming from the truncated version of the 
Hamitonian, as explained in the next section. 

3. Beyond mean field 

The first step in going beyond the mean field description is the introduction of the 
displaced boson operators, as in [T3J |3] , 

c ] = b ] - p c , c = b- p c , (7) 

which are defined to satisfy the property c\a c ) = 0. While in the normal phase, where 
p c = 0, there is not displacement, using c' , c allows for a general treatment of the Lipkin 
model beyond mean field. 

If the ground state is well described by the coherent state (the vacuum in the normal 
phase), it is valid to make the approximation 

It must be stressed that this approximation fails badly close to the phase transition, 
where the truncated Hamiltonian cannot be employed. Far from the phase transition 
region it is possible to expand the square roots in Eq. ^ in series, conserving terms 
of order N 1 , N2 and N°, and neglecting all negative powers of N. Care must be taken 
because p c is of order N. The truncated version of the Hamiltonian iJ, which is 
quadratic in the new bosons, is 

H® = A + Bc^c + C {c^ 2 + c 2 ) (9) 
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A = -7c A 



+ 



-7 C + 



lx 

AN 

lx 



(N - 3p 2 c + Ap\ (N 



pI 



+ 



7^ 
AN 



{n-p: 



(10) 



with 

2N 2N 
(N-5pp lx {N-pp ly 

AN AN 

There is an extra term proportional to c* + c whose coefficient vanishes exactly when 
employing the values of p c given in Eq. The Hamiltonian can be diagonalized 
through the Bogoliubov transformation (TO] 



C 



,t 



cosh 



-©- 

.2". 


+ sinh 


e- 

.2~_ 


a, c = cosh 


e- 

.2". 


a + sinh 


-0" 
.2". 



in terms of the new bosons a\ a which, when replaced in Eq. ([9]), it reads 

-0,2 



H (t) 



+ Csinh[e] + (Bcosh[6] + 2Csinh[0])a t a 



A + B sinh 

■ sinh[6] + Ccosh[eA (a t2 + a 



+ 



(12) 



The last term cancels out by selecting tanh[0] 
the final diagonal form [10] 



-2^. The truncated Hamitonian has 



H {t) = A+ - (^B 2 -AC 2 - l) + VB 2 -AC 2 a j a = Ne$ + Aa j a. (13) 

The ground state of the truncated Hamiltonian is the new vacuum |0), which satisfies 
a|0) = 0. The coefficient of the last term is the gap A, the energy separation between 
the ground and the first excited state in the normal phase, and between the ground and 
the second excited state in the deformed region, where the first excited state becomes 
degenerate with the ground state in the thermodynamic limit [TUHIS]- Its explicit form 
is 

\j{lx -7c)(-7c + 7y) region I , 

region II . 



A = VB 2 -AC 2 



(14) 



{lx -7c) (lx-ly)/lx 

Figure [2] plots the g function of ^ x as dashed blue line. It becomes null at the 

phase transition r y x = j c . The excitation energies of the first and second excited states, 
obtained through exact diagonalization for N=40, are also displayed. The truncated 
Hamiltonian allows a good description of the gap, which becomes exact when N — > oo. 
The minimum of the excitation energy is a precursor of the phase transition at finite N, 
which takes place at a different value of 7 X for each iV [22] . 

The ground state energy per atom i s given by the constant terms in 
Eq. (ph, divided by N 



r(*) 
"9s 



7c I 7e I 

2 2N 2N 



' 2N 



jj\j{lx-lc){l y -lc) region I , 

+ 2N^f~x - 7c) (lx - ly)/lx region II. 



(15) 
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Figure 2. (Color online) The gap, Eq. (14), as function of -f x , shown as a dashed 
blue line for 7 C = — 1 and j v = 1, and the first and second excitation energies, obtained 
numerically for N= 40, displayed as continuos green lines. 



It is plotted as a continuous blue line in Fig. [T] (left), and displays a spike at y x = j c 
which is spurious, an artifact of the truncation, which vanishes as N — > oo. 

The fraction of excited atoms is obtained by expressing J z in terms of the new 
bosons {a\ a), and evaluating n e = "^4r- + 1 in the new vacuum. 



N 

7a:+7«— 2 7c J_ 



region I , 



n(t) , 2i V^-7c)(7,-7c) " 

' + I ;v T - h region II. 

The curve n$ vs. 7 X is displayed as a continuous blue line in Fig. [I] (right). It follows 
closely the mean field prediction, except in a vicinity around 7 X . = j c , where it diverges. 
In this case the pathologic behaviour survives in the thermodynamic limit. It reminds 
us that in this region the solutions of the truncated Hamiltonian are no longer valid. 
There are, however, reliable ways to extract useful information from the spurious results 
obtained around the phase transition. They are the subject of the next section. 



4. renormalisation and critical exponents 

The renormalisation procedure postulates that it is possible to extract the correct 
functional form 0(7, N), smooth and finite for any finite N, of any observable ( -*- ) (7, N) 
which becomes singular at the phase transition [23] when described using the truncated 
Hamiltonian. IfO (t) (7,iV) is singular as N^{^j — 7 C ) a , employing a scaling function 
[A^ y (7 — 7c)] Q the regular function is built as [10J 

0(j,N) = OW( 1 ,N) [AP(7- 7c )] Q 

iV^(7 - i c y a [JV( 7 - lc)T = (17) 

where the power a in the second term was selected to cancel out the singularity at y c , 
and v defines the class of universality to which the model belongs. In the Lipkin model 
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numerical analysis points to v — | (6J E21 ES] . We will show that this number can be 
deduced analytically using the fidelity susceptibility. Also we will exhibit the singular 
behaviour of the point ( 7a; , j y ) = ( 7c , 7c ) in the parameter region, where v has a different 
value when approached along the lines j y = 7 C or j x = 7c . 

The singular term in the energy per atom, which gives rise to the spurious spike 



when 7-e — > 7 C , behaves like a/7x — 7c/ 'N, as shown in Eq. (15). Renormalizing it we 
obtain 

1 



^r(7„iV)->(7,-7c) 5 iV- i AT3( 7x -7 c ) 



4 



Ns. 



The singular term can be calculated numerically, substracting from the exact energy 



per atom the regular terms in Eq. (15). It goes to zero as function of N as predicted 
The divergence in the fraction of excited atoms can be manipulated in the same 



way. Taking the divergent term from Eq. (16), it gives 



2 



< en (7x,iV)^(7,-7c)~ 5 A r " 1 [iV3( 7:E - 7c )j 2 =AT-3. (19) 
In Ref. [TU] this dependence on N was also confirmed numerically. 

4-1- Fidelity susceptibility 

In classical information theory the fidelity measures the accuracy of a transmission [24J . 
It also provides a powerful tool to study quantum phase transitions [2S]- For a pure 
state which varies as a function of a control parameter A, the fidelity is defined 

as 

F(\,\ + 5\) = m\)m\ + 5\))\ 2 . (20) 

It diminishes when there are sudden changes in the wave function, which for finite N 
represent a precursor of a phase transition, and has a minimum where these changes 
are the largest. Even more sensitive is the fidelity susceptibility [26], closely related to 
the second derivative of the fidelity. It has a maximum reflecting the phase transition. 
Expressing the Hamiltonian under study as H = Hq + Aifj, which makes explicit the 
control parameter, the fidelity susceptibility can be calculated as [251 



(21) 

where |^(A)) denotes the k eigenstate of H with energy £fc(A), and k = refers to the 
ground state. 

In our case, we select Hj = J%/N and study the dependence on j x . The way to 
proceed is to express it in terms of the bosons (V , a, as 

Jx = Ji + ha 1 a + j 3 (a* + j 4 (a t2 + a 2 ) . (22) 
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The first two terms do not contribute to xf- The third one connects the ground state 
with a one boson state, with energy E\ = A, the last one with a two boson state with 
energy E 2 = 2A. In the normal phases, region I, we find 



* = °- * = £ vK (23) 



and in the deformed phase, region II, 



( 7 2_ 7 2)( 7x _ 72/ ) 



47x V 7* 



;d _ iV ( 7c 2 (5 TsB - 37„) - 7c7. (T« + 37,0 + 27*7*) 



%lx\jlx{ll -7c) S -7s,) 



Substituting in Eq. (21) we obtain 




Xl 32 ( 7x - lc f ' X ^ 4 7, 3 ^ (7c + 7 J(7, - 7.) ' (26) 

Note that in the deformed region II 7a . < 7c = — e < and 7x < 7j) , making the argument 
in the square roots, (j|) 2 and Xf positive. 

The fidelity susceptibility xf exhibits very interesting features. It is divergent in 
both phases as 7x — > 7c , but with different powers of ( 7x — 7c ) and of N in each phase 
Renormalizing the fidelity susceptibility the following exponents are found 





a 


P 


/3 + a v 




>7c 


2 





2v 


Ix 


< 7c 


i 

2 


1 


1 + - 

x ~ 2 



It implies that the renormalized fidelity susceptibility xf scales as N 2v in the 
normal region, and as iV 1+ ^ in the deformed region. As xf is continuos at the 
phase transition for any finite N, the two exponents should be equal, implying that 



" = !■ 



In this way we have shown analytically that the Lipkin Hamiltonian belongs to this 

4 

class of universality. It follows that xf diverges as iV3 at the phase transition. This 
dependence of xf at 7c is confirmed numerically [22J. It is shown in the left of Fig. [3j 
where the maximum value of xf is plotted against iV in a log 2 -log 2 plot. The points are 
clearly along a straight line, with slope 1.37 ~ 4/3. The numerical fit is xf — 0.175 iV 1,37 . 
The values of N range from 2 10 = 1024 to 2 16 = 65536. 

The above deduction is valid for j y > 7c in the normal phase, region I, and for 
ly < lx in the deformed phase, region II. There is a very singular behaviour along the 
line jy = 7c in the deformed region II. In this case the fidelity susceptibility has the 
functional form 

rtb* = ^ = -r^^i^x^^- (27) 

(7* - 7c) 47 3 V 7x + 7c 
Numerical studies [22] have shown that u(j y = j c ) = 1, and xf — > N 2 as it approaches 
the phase transition. This behaviour is plotted in the right of Fig[3j with slope 2.001. 
The numerical fit is xf = 0.892 iV 2 along the line 7 y = 7c , ■y x — > 7 C _. 
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Log,N 



Figure 3. (Color online) Log-log plots of the maximum of the fidelity susceptibility 
as function of the number of particles N in the system, for j c = — 1 and j y = 1.0 (left) 
and 7 y = —1.0 (right). 



5. Conclusions 



We have studied the Lipkin model of collective spins employing the Holstein-Primakkof 
mapping of the quantum rotor to a bosonic field. We have shown how the different 
phases of the system are found by using the coherent Heisenberg-Weyl state, which is 
an eigenstate of this boson. The possibility of going beyond the mean field description 
was reviewed. It was obtained by introducing a second set of displaced bosonic operators 
whose vacuum is the coherent state, and expanding the square roots appearing in the 
mapping of the angular momentum operators as an infinite series, replacing it in the 
Hamiltonian, and truncating by keeping terms of order N 1 , N 1 ^ 2 and N°. In this way 
a truncated Hamiltonian is built, which is quadratic in the new bosons, and can be 
diagonalized exactly by means of a Bogoliubov transformation. It provides a harmonic 
description of the Hamiltonian, with an excitation energy, the gap, which vanishes at 
the phase transition. It was shown that this description is consistent with the exact 
behaviour of the system at finite number of spins N, obtained through a numerical 
diagonalization of the full Hamiltonian. 

On the other hand, it was also exhibited that the predictions for the ground state 
energy per particle and for the fraction of excited atoms exhibit a spike around the phase 
transition which is spurious, an artifact of the truncation, which is itself not longer valid 
in a vicinity of the phase transition. The size of this spike goes, in the case of the energy 
per particle, to zero in the thermodynamic limit, while for the fraction of excited atoms 
is infinite, and remains infinite even in the limit N — > oo. These singularities, which 
inhibit the use of the truncated Hamiltonian close to the phase transition, carry anyway 
useful information. It was obtained through renormalisation, a recipe to built smooth 
and well behaved functions from singular ones by multiplying them with an appropriate 
function of the control parameter and the number of particles. 

The behaviour for large N of the singular term in the energy per particle and of the 
number of excited atoms was obtained in this way. It was also shown that the fidelity 
susceptibility, an observable widely used in quantum optics, provides not only a very 
efficient tool to detect the precursor of the phase transition at finite N, but allows one 
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to obtain analytically the class of universality to which the Lipkin model belongs: the 
one which is renormalized with a general factor (iV 2 / 3 (7 — 7 c ) a . The exponent 2/3 was 
obtained asking for the fidelity susceptibility to be continuos at the phase transition, and 
taking advantage of the fact that it has different functional forms in the different phases 
of the system. It was also show that there is a specific region, the line where 7 y = 7 C 
and 7^ < j c , in which the renormalizing function is different, i.e. the system behaves in 
a different way along this line. Numerical calculations confirm the N dependence of the 
fidelity susceptibility in both regions. 

We hope that this review of the subject would serve both as an introduction to the 
powerful techniques and to the care that must be taken in its use. 

We thank J. Vidal for his comments. J.G.H. and E.N-A. thank P. Domokos and D. 
Nagy for their hospitality and valuable conversations. This work was partially supported 
by CONACyT-Mexico and DGAPA-UNAM (projects IN102811 and 103212). 
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